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Abstract. The optical spectrum of the infrared source 
IRAS 04296+3429 (optical counterpart - GO la star, V = 14.2) 
was obtained with the echelle spectrometer PFES at the prime 
focus of the 6 m telescope. We discover emission bands (0,0) 
and (0,1) of the Swan system of the C2 molecule in the optical 
spectrum of IRAS 04296+3429. Comparison with the spectrum 
of the Hale-Bopp comet leads us to propose that in both cases 
the same mechanism (resonance fluorescence) is responsible 
for the emission in the C2 molecular bands. 

Several strong absorption features whose positions coin- 
cide with known diffuse interstellar bands (DIBs) are revealed 
in the spectrum of IRAS 04296+3429. 

The infrared spectrum of IRAS 04296+3429 shows the fa- 
mous 21 fim feature (Kwok et al. 1989), but this object has 
not been observed by KAO (Omont et al. 1995). However, like 
IRAS 051 13+1347, IRAS 05341+0852 and IRAS 22223+4327 
(Kwok et al. 1995, Szczerba et al. 1996), our detailed mod- 
elling of its spectral energy distribution suggested that this 
source also should show the 30 fim band. In fact, ISO dis- 
covered a broad, relatively strong feature around 30 fim for 
IRAS 04296+3429 (Szczerba et al. 1999). 

The surface chemical composition of the source 
IRAS 04296+3429 is metal-deficient (the averaged value 
of the abundances of the iron group elements Ti, V, Cr and Fe 
relative to the solar values is [M/H] Q = —0.9) and has been 
considerably altered during the evolution: carbon, nitrogen and 
s-process elements are overabundant relative to the metallicity. 

The totality of physical and chemical parameters derived 
for IRAS 04296+3429 confirms a relation between presence of 
the feature at 21 /im in the spectrum of a carbon rich star and 
an excess of the s-process elements. 

Key words: Stars: abundances - Stars: atmospheres - Stars: 
AGB and post-AGB - Stars: individual: IRAS 04296+3429 



1. Introduction 

This paper is a continuation of our series of papers on 
the investigation of stars which are believed to be in the 

Send offprint requests to: V.G. Klochkova (valenta@sao.ru) 



post-Asymptotic Giant Branch (post-AGB) stage of evolu- 
tion (Klochkova 1995; Zacs et al. 1995, 1996; Klochkova & 
Panchuk 1996a, 1998; Klochkova et al. 1997a; Klochkova & 
Mishenina 1998). The post-AGB stars (hereafter also referred 
to as proto-planetary nebulae - PPNe) being in the transi- 
tion phase from AGB to planetary nebulae offer an opportu- 
nity to study in detail a chemical composition which has un- 
dergone changes due to nucleosynthesis and mixing processes 
in the course of the stars evolution. Here we present new re- 
sults for the peculiar supergiant with a large infrared excess 
IRAS 04296+3429. 

On the 12/25/60 /im colour-colour diagram from the 
IRAS data infrared source IRAS 04296+3429 (hereafter 
IRAS 04296), associated with a faint carbon-rich star (Omont 
et al. 1993; Loup et al. 1993) classified as type GO la by Hriv- 
nak et al. (1994), is located in the region occupied by planetary 
nebulae, non-variable OH/IR stars, and proto-planetary nebu- 
lae (Iyengar & Parthasarathy 1997). The object IRAS 04296 
belongs to the small group of sources which show a spec- 
tral feature around 21 /im (Kwok et al. 1989). This feature is 
seen only for some post-AGB objects and has not been de- 
tected either in the preceding (AGB) nor in succeeding (PN) 
evolution stage. Note that a search for new 21 /im emitters by 
means of ISO SWS observations among candidates selected by 
Henning et al. (1996) failed to give any detections (Henning, 
private communication). Using a medium-resolution (3 A) op- 
tical spectrum, Hrivnak (1995) found that IRAS 04296 is a 
strongly reddenned (E(B-V)= 1.3) G-star with features indica- 
tive of high luminosity, with molecular absorption features of 
C2 and rarely observed features of C3 of c ire urns te liar ori- 
gin and quite strong absorption lines of s-process elements 
(Ba, Sr, Y) indicating that outer layers of the atmosphere of 
IRAS 04296 have been enriched by products of nucleosynthe- 
sis. Therefore this star is very well suited for the study of de- 
tailed chemical abundances which have been changed by the 
third dredge-up. Indeed, Decin et al. (1998), using high res- 
olution spectra, obtained the chemical abundance pattern for 
this object and concluded that its metal-deficient, carbon-rich 
atmosphere has large overabundances of s-process elements. 

In Sect. 2 we describe our observational material for 
IRAS 04296 and discuss its molecular features, comparing 
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them with the corresponding spectrum of the Hale-Bopp 
comet. Sect. 3 is devoted to presentation of the main param- 
eters and detailed analysis of the chemical composition of 
IRAS 04296 derived from our optical spectra. The next sec- 
tion presents modelling of spectral energy distribution for this 
source with the aim to get insight into its physical parameters 
(mainly to determine the stellar effective temperature which is 
crucial for the chemical composition estimation). Finally, in 
Sect. 5 we discuss the results obtained and compare them to 
the results for related objects. 

2. Optical spectrum of IRAS 04296 

2.1. Observations and spectra reduction 

We have obtained spectra of IRAS 04296 with the CCD 
(1140x 1170 pixels) equipped echelle spectrometer PFES 
mounted at the prime focus of the 6 m telescope of SAO RAS 
(Panchuk et al. 1998). The echelle-grating with 75 gr/mm and 
with a blaze angle of O = 64.3 ° was used. A diffracton grat- 
ing with 300 gr/mm was used as the cross-disperser. The cam- 
era has f = 140 mm. The projected angular size of the input slit 
is 0.54 arc sec. 

We observed IRAS 04296 on October 07, 1996 
(JD2450363.6) and February 26, 1997 (JD2450506.3). 
The echelle frames with 25 echelle-orders cover the spectral 
region AA 4420-8300 A. The average spectral resolution was 
0.4 A. The signal-to-noise ratio was in the range 50-110 for 
different spectral orders. 

All usual procedures needed for echelle-images reduction 
(bias subtraction, cosmic ray removal, optimal order extraction, 
rebinning) were made using the ECHELLE context of the MI- 
DAS system. A Gaussian function approximation was made for 
the measurement of equivalent widths. The comparison spec- 
trum source was an argon-filled thorium hollow-cathode lamp. 

The distinctive features of the optical spectrum of the 
source IRAS 04296 are a peculiar profile of the H Q line (see 
Fig. 0), molecular emission bands and very strong absorption 
lines of ionized atoms of s-process elements (Y, Zr, Ba, La, 
Ce, Pr, Nd). For example, the equivalent widths of Ball lines 
(6141 and 6496 A) exceed 0.6 A. 

2.2. Emission molecular bands 

Absorptional bands of several molecules (C2, CN, TiO, etc.) 
are often present in the spectra of post-AGB stars (see, for ex- 
ample, Hrivnak 1995; Bakker et al. 1997). However, molecular 
emission features are only very rarely observed in the optical 
spectra of PPNe. One such example is RAFGL2688 (the Egg 
Nebula) for which Crampton et al. (1975) observed emission 
features of the C2 molecule in a medium resolution spectrum. 
On the other hand, it is well known that cometary nuclei spectra 
show prominent Swan band emission. 

In both spectra of IRAS 04296 we have discovered strong 
emission in the (0;0) and (0;1) bands of the Swan system of 
the C2 molecule. On Figs. we present a comparison be- 
tween the spectrum of IRAS 04296 (observed on February, 26, 
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Fig. 1. The IRAS 04296+3429 spectrum near H«. 



1997) and that of the Hale-Bopp comet (observed on March 
30, 1997 with the same spectrometer) around bands (0;1), (0;0) 
and (1;0), respectively. From Figs. ||-0it is clear that emission 
band (1;0) at 4735 A is absent in the spectrum of IRAS 04296 
while the bands (0;1) at 5635 A and (0;0) at 5 165 A are re- 
liably measured. Hrivnak (1995) obtained the spectrum of 
IRAS 04296 inside the blue spectral region, 3872^1870 AA, 
therefore he could not observe emission features of C2 at 5165 
and 5635 AA. 

To understand the observed ratios between different bands, 
we have estimated the temperature function for monochromatic 
coefficient of absorption per molecule (a\) for the Swan bands 
in the "just overlapping" approximation (JOA, Golden 1967). 
For the microturbulent velocity £ t = 7 km/sec this approxi- 
mation works well near the band heads. Values of a\ for band 
heads of (0;0) at 5165 A and (1;0) at 4735 A coincide within 
0.2 dex for the temperatures range 3000-7000 K, while o\ for 
the band (0;1) at 5635 A is systemically lower by about 0.6 dex. 
Taking into account these relations between a\& for different 
band heads and since we do not observe the 4735 A band in 
RAFGL2688 and IRAS 04296, we can conclude that it is im- 
possible to describe the intensity ratios of C2 emission bands 
for these objects by means of an equilibrium vibrational tem- 
perature in the 3000-7000 K range. 
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Fig. 2. Comparison of the IRAS 04296 + 3429 C 2 Swan band 
head (0; 1) at A = 5635 A with that for the Hale-Bopp comet nu- 
clei 



Fig. 3. The same as Fig. 2 but for C2 Swan band head (0;0) at 
A = 5165A 



To explain emission bands intensities for comets the 
mechanism of resonance fluorescence has been proposed 
(Zanstra 1928, Swings 1941). In that case population of 
vibration-rotational levels for the molecule is described by the 
Boltzmann approximation, however the value of T in the expo- 
nent no longer has the meaning of equilibrium temperature but 
it is a distribution parameter only. We suggest that the same 
mechanism could be responsible for the observed emission 
bands of IRAS 04296. However, it is clear from Figs. ||-[|that 
there are significant differences in the equivalent widths of the 
emission bands in the spectra of our supergiant and of the Hale- 
Bopp comet nucleus. They could be explained by a difference 
of radiation fluxes which illuminate C2 molecules in these ob- 
jects. The temperature of IRAS 04296 (T eff around 6300 K) 
is sufficiently higher than that for the Sun, therefore the band 
(1,0) at 4735 A for IRAS 04296 should be stronger than that for 
the Hale-Bopp comet nuclei. However, Fig.^ shows the oppo- 
site behaviour. It could mean that radiation field of IRAS 04296 
which excites the C2 molecules is strongly reddened by matter 
located between its photosphere and the region which produces 
the C2 emission. 

Together with the emission bands of the Swan system 
(Klochkova et al. 1997b) absorption bands of the Phillips sys- 
tem (1:0), (2;0), (3;0) have been revealed in the spectrum of 



IRAS 04296 (Bakker et al. 1997). Let us try to explain this 
phenomenon within the resonance fluorescence mechanism or- 
dinary used to interpret comets' spectra. As a first approxima- 
tion, we assume that the vibrational distribution corresponds to 
the effective temperature of the star illuminating a circumstel- 
lar envelope if vibrational transitions in the low triplet state of 
a homonuclear molecule are strictly forbidden. But even when 
interpreting comets' spectra such an approach appears to be 
too poor. The intensity distributions for different systems of 
bands and for bands of individual systems of the resonance 
fluorescence of the C2 molecule have been considered in pa- 
pers by Krishna Swamy & O'Dell (1977, 1979, 1981). The 
intensities of bands have been calculated taking into account 
the excitation of the Swan, Ballick-Ramsay and Fox-Herzberg 
triplet systems, Phillips and Milliken singlet systems as well 
as singlet- triplet transitions in low states. It has been shown, in 
particular, that at the value of the moment of singlet-triplet tran- 
sitions |R C | 2 = 10~ 5 and at the heliocentric distance of a comet 
d=l a.u. the ratios of intensities of sequences At/ = 0, 1, —1 
in the Phillips system to the intensity of sequence Av — of 
the Swan system is equal to 0.094, 0.11 and 0.04, correspond- 
ingly (Krishna Swamy & O'Dell 1981). This agrees well with 
results of measurement of comets' spectra. Using these results 
of Krishna Swamy & O'Dell (1981), we may suppose that the 
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Fig. 4. The same as Fig. 2 but for C2 Swan band head (1;0) at 
A = 4735A 



intensity of main bands of the Swan system is ten times higher 
than that in the Phillips system. 

Now consider the case of IRAS 04296. Let us add such 
an emission spectrum of the C2 on the stellar continuum. In 
order to observe the emission bands of both the Swan and 
the Phillips systems over the continuum in such a combined 
spectrum, the stellar flux at A = 5 165 A must be at least 10 
times higher than near A = 7720A. From Kurucz's (1979) ta- 
bles it follows that the ratio of the fluxes near these wave- 
lenghts for the Sun (the emitter in the case of comets) is equal 
to F A 5i7 5 /F A 775o = 1.5. For the model with T cff = 6300 K this 
ratio is equal to FA5175/FA7750 = 1-9. From the real spectral 
energy distrubution observed for IRAS 04296 (Kwok 1993) the 
ratio of the fluxes is essentially smaller: FA5165/FA7720 = 0.1. 
Therefore, the conditions to observe the absorption bands of 
the Phillips system and the emission bands of the Swan system 
may arise inside the circumstellar envelope of IRAS 04296. 

2.3. Radial velocity of the IRAS 04296 

We measured the radial velocity V r using our best 
spectrum, that obtained in February, 1997. The average 
value of the radial velocity from numerous metal lines 
(V r = —56.0 ± 0.8 km/s) and from the H Q absorption line 
(V r = — 54.5km/s) agree within the accuracy of measure- 



ment. This radial velocity is consistent with a membership 
of an old population as suggested by the low metallicity (see 
Table 1). It should be noted also that the value of radial ve- 
locity we derived agrees with the value of V r = — 59km/s 
which was given by Omont et al. (1993) from CO data and 
V r = -62 km/s ±1.0 km/s by Decin et al. (1998) from op- 
tical spectrum. There is still no sign of an essential temporal 
variability of V r for the object. 

2.4. Absorption bands identified with DIBs 

From comparison of observed and synthetic spectra of 
IRAS 04296 we discovered some strong absorptional features 
whose positions coincide with known diffuse interstellar bands 
(DIBs) (Jenniskens et al. 1994). In Fig.|| we illustrate the pres- 
ence of DIB's by showing a spectral region of IRAS 04296. We 
have calculated the synthetic spectrum using the code STARSP 
(Tsymbal 1995) and the atmospheric parameters and abun- 
dances of chemical elements we here obtained. It should be 
noted that for such a comparison in the spectral range near A 
6270-63 10 A the telluric spectrum has been removed from the 
observed spectrum. 

In the following paper we plan to study in detail such iden- 
tified with DIB's absorptions, we have revealed in the spec- 
tra of several related objects (IRAS 04296, IRAS 23304+6147, 
IRAS222223+4327), here we limited ourself by such short in- 
formation. 



3. Determination of atmospheric parameters and 
calculation of chemical composition 

For understanding of an object at an advanced evolutionary 
stage, it is very important to know its metallicity and detailed 
chemical abundance pattern. Our echelle spectra provide such 
a possibility due to their large wavelengths coverage. 

To study the chemical composition, we have used the plane- 
parallel homogeneous models generated by the MARCS pro- 
gram (Gustafsson et al. 1975). It should be noted, however, 
that unstable and very extended atmospheres of supergiants 
probably require more advanced model atmospheres. There- 
fore, our results should be treated as only preliminary ones. 
For a chemical composition calculation by the model atmo- 
sphere method, one needs to know the values of the effec- 
tive temperature (T e ff), surface gravity (logg) and microtur- 
bulent velocity (£ t ). Determination of T c ff is problematic even 
for normal supergiants due to their extended atmospheres and 
significant non-LTE effects. In the case of so peculiar a su- 
pergiant as IRAS 04296, for which the energy distribution is 
strongly distorted by interstellar and circumstellar extinction, 
determination of T e ff is the most difficult problem. We can- 
not use for this purpose equivalent widths and profiles of H I 
lines (well known criteria of atmospheric conditions for nor- 
mal supergiants), since these lines are strongly distorted in the 
spectrum of IRAS 04296 as seen in Fig . [l]. 

Therefore, we have applied the spectroscopic method for 
temperature determination of IRAS 04296, forcing the abun- 
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dance derived for each line to be independent of the lower ex- 
citation potential (EP). We have estimated that T e g = 6300 K 
with an internal uncertainty A T c g = 250 K. To check the re- 
aliability of our determination we have modelled the spec- 
tral energy distribution for this source (see Sect. 4) and got 
a very similar temperature near 6500 K. The surface grav- 
ity log g = 0.0 was estimated through the ionization balance 
of the Fel and Fell abundances. The errors on the parame- 
ter log g is determined by forcing a maximum difference be- 
tween e (Fe I) and e (Fe II) to be 0. 1 dex (where here and here- 
after, loge(X) = logN(X) -logN(H)). It should be noted 
that the hydrogen abundance logN(H) = 12. Such a difference 
is achieved by varying the logg value by ±0.2 keeping other 
parameters (T e fi and £t) constant. The microturbulent veloc- 
ity value based on equivalent widths (W) of Fel and Fell 
lines is quite high, equal to 7 km/s. This value is determined 
with an uncertainty of ±1.0 km/s, which is typical for F, G- 
supergiants. 

To illustrate the choice of model parameters for the ob- 
ject IRAS 04296 in the Fig.|| are shown the excitation poten- 
tial - abundance diagram and the equivalent width - abun- 
dance diagram for lines of neutral (dots) and ionized (crosses) 
iron atoms. As follows from this figure, there are no essen- 



tial dependences for values considered. The large dispersion 
is mainly explained by errors of measurement of equivalent 
widths of weak absorption lines for such a faint object as the 
IRAS 04296 (see, for example, the similar dispersion on the 
Fig. 1 in the paper by Decin et al. (1998) for the brighter object 
IRAS 22223+4327, V=9.7). 

We have checked the determination of IRAS 04296 model 
parameters using weaker Fel and Fell lines and concluded that 
the parameters are steady within the erorr box up to W = 100- 
150 mA. This can also be seen from Fig.^. 

It is well known that the plane-parallel static model at- 
mosphere method does not give correct abundances for high 
luminosity stars (luminosity classes la, Ia+). The profiles of 
the spectral lines observed are broadened by non-thermal 
mechanisms whose influence may be variable at different 
levels in the atmosphere. Therefore, to obtain more reli- 
able estimates of chemical element abundances we use weak 
lines with W <250mA. The average values of the equiva- 
lent widths W we used for the abundances calculations are 
also given in Table 1. Only the Ball abundance was calcu- 
lated using 3 very strong lines: W(A 5853.67) = 464mA, W(A 
6141.71) = 679 mA and W(A 6496.90) = 738 mA, because the 
weaker lines of this element were not available. In general, the 
weak lines formed in deeper atmospheric layers are more cor- 
rectly described by the standard static model. The limitation 
of equivalent width of lines used to W <250mA significantly 
reduces the influence of uncertainty in the choice of £ t - Note, 
however, that the main factor in the abundance errors for most 
species remains the uncertainty of the T c ft value. Therefore, 
we have checked our estimation of T ff by modelling of spec- 
tral energy distribution for IRAS 04296. 

Computed abundances of 26 chemical elements are pre- 
sented in Table 1 . In the head of the Table 1 parameters of 
the adopted model atmosphere are shown. The dependence 
of chemical composition determination on uncertanties of the 
model atmosphere parameters is discussed in Zacs et al. (1995). 
In the second column of Table 1 derived abundances are given 
as loge (X), while in the third column estimated uncertainties 
of a = A log e (X) are shown. In the next column, the number 
of spectral lines used for chemical composition calculation is 
indicated. 

A lot of absorption lines of different elements (CNO- 
elements, light metals, iron group elements, Ce, Nd, Eu) have 
been reliably measured in the spectrum of IRAS 04296. It is 
important that we have not found any dependence of the abun- 
dances of these species on the equivalent width or on the exci- 
tation potential. Therefore the microturbulent velocity does not 
vary between different chemical elements. 

The gf-values for most of the spectral lines used for 
the abundance calculations were taken from the list used by 
Luck (1991). The S and CNO-abundances were determined by 
using the gf-data from Waelkens et al. (1991) and Giridhar et 
al. (1994). The list of lines with the adopted gf-values, exci- 
tation potentials of the lower level and equivalent widths we 
measured for the object IRAS 04296 are available by e-mail 
(valenta@alba.sao.ru). 
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Table 1. Model atmosphere parameters adopted and abundances of chemical elements. Here, n is number of lines used for 
calculation, a - the standard deviation, W - the average equivalent width, in mA, of lines used for the content calculation 
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To verify the method of analysis we observed with the same 
spectral device the normal supergiant a Per. The same pro- 
cedures for processing and the same list of lines were used 
for analysis of the a Per spectrum. This supergiant, whose 
parameters, T e s = 6500 K, logg = 1.2, £ t = 4.7 km/sec are 
very close to the object studied, is very convenient as a stan- 
dard for the method testing because of its membership in the 
young open cluster aPer which has solar chemical composi- 
tion (Klochkova & Panchuk 1985; Boesgaard 1989). Using its 
membership of this cluster, we may predict that a Per also has 
normal solar chemical composition (aside from the expected 
nonsolar CNO triad abundances relative to iron). As it is shown 
in Table 1 aPer has indeed the abundances of chemical ele- 
ments close to solar ones, except for CNO and several elements 
whose abundances are calculated with a large uncertainity due 
to a small number of spectral lines used. 



4. Spectral energy distribution 

Details of the computer code used for solution of the radia- 
tive transfer in dusty envelopes can be found in Szczerba et 
al. (1997). In brief: the frequency-dependent radiative transfer 
equation is solved for a dust under assumption of spherically- 
symmetric geometry for its distribution taking into account par- 
ticle size distribution and quantum heating effects for the very 
small dust particles. 

The modelled source is certainly C-rich (see Omont et 
al. 1995). Therefore, for modelling of its spectral energy dis- 
tribution (SED) we assumed that dust is composed of: poli- 
cyclic aromatic hydrocarbons (PAH) for dust sizes a between 
5 and 10 A (see Szczerba et al. (1997) for details concerning 
PAH properties), amorphous carbon grains (of AC type from 
Rouleau & Martin 1991) for a > 50 A, and dust with an opac- 
ity obtained from averaging of the absorption efficiences for 
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Fig. 6. Upper: iron abundance Fel (circles) and Fell 
(crosses) calculated for IRAS 04296 with model parameters 
T e ff = 6300 K, logg = 0.0 and £ t = 7.0 km/sec using lines 
with different EP of a low level; bottom: the same as a func- 
tion of the equivalent width W 



PAH and AC grains according to the formula: 
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)PAH 

' abs, 
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for grain sizes between 10 and 50 A. Here: f = 1 for a= 10 A 
and f = for a = 50 A. Dust with opacity values constructed in 
this way allow us to use a continous distribution of dust grain 
sizes and fill the gap between properties of carbon-bearing 
molecules and small carbon grains. 

The 21 pm feature was approximated by a gaussian with 
parameters determined from modelling of IRAS 07134+1005 
(centre wavelength equal to 20.6 /im, and width of 1.5 pm) 
which has the strongest feature among the known 21 /im 
sources. In the case of 30 pm band we used the addition of two 
half-gausians with the same strength and different width. Ini- 
tial fit was done to IRAS 22272+5435 and its parameters were: 
width for short wavelength side ctl = 4 /im, width for long 
wavelength side ctr = 9 /im and central wavelength 27.2 pm 
(see Szczerba et al. 1997). For modelling of IRAS 04296 we 
have reduced the strength of this feature by 50%. Superposi- 
tion of the 21 and 30 pm features was added to the absorption 
properties of amorphous carbon in order to construct an empir- 
ical opacity function (EOF). 

In Fig. ^ the best fit obtained from the solution of the radia- 
tive transfer problem including quantum heating effects for the 
PAH grains is shown together with observational data which 
will be described in detail elsewhere. Note, however, that we 
present also two sets of photometry (from B to M band) cor- 
rected for interstellar extinction (open symbols) according to 
the average extinction law of Cardelli et al. (1989), assuming 
that total extinction at V is 1.0 or 2.0 magnitudes and plot- 
ting only the smallest and largest value of corrected fluxes 
at given band. This estimate of the total extinction range can 
be inferred from the analysis of data presented by Burstein & 
Heiles (1982). 

The best fit to the spectral energy distribution of 
IRAS 04296 is shown by heavy solid line (see Table 2 for de- 
tails concerning parameters of the model). Our modelling pro- 
cedure was such that we tried to get fits to the SED which 
fall in between the extinction corrected fluxes. In this way, we 
have taken into account not only the effect of the circumstel- 
lar extinction but also of interstellar extinction. The thin long- 
dashed line represents the input energy distribution of the cen- 
tral star for log g = 0.5 and T c ff = 6500 K according to model at- 
mosphere calculations of Kurucz (private communication). The 
heavy short-dashed line shows the fit which was obtained with 
the same assumptions but changing the effective temperature of 
the star to 6000 K. As one can immediately see in the IR range 
of the spectrum the quality of the fits are very similar. However, 
in the optical and ultraviolet (UV) part of the spectrum the fit 
assuming T c ff = 6000 K is not able to explain extinction cor- 
rected data. In consequence, we are quite convinced that our 
estimation of T off for IRAS 04296 close to 6500 K is reason- 
able and, what is even more important, agrees pretty well with 
the spectroscopic estimation (6300 K). Note that spectral type 
of this source was found to be GO la from the low resolution 
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Fig. 7. Fit to the spectral energy distribution of 
IRAS 04296 + 3429 obtained with empirical opacity func- 
tion, taking into account quantum heating effects (heavy solid 
line) and assuming that star radiates according to a model 
stellar atmosphere calculations for log g = 0.5 and T c ff = 6500 
K (thin long-dashed line). The heavy dashed line shows the 
fit obtained for star with T c h =6000 K. The thin solid line 
underlying 21 and 30 /im features represents the estimated 
continuum level 

spectrum (Hrivnak 1995) which implies an effective tempera- 
ture of around 5500 K for the star if we asume that the same 
relationship applies for post-AGB supergiants as for "normal" 
ones (see Schmidt-Kaler 1982). For such a low temperature we 
were not able to fit even the reddenned data in the UV. 

The thin solid line in the wavelength range from about 1 8 
to 48 fim represents the model continuum level found after so- 
lution of radiative transfer equation for dust without EOF using 
the parameters as in Tab. 2 while keeping the dust temperature 
(or probability distribution of dust temperature) the same as for 
the case of dust with EOF. Taking into account the estimated 
continuum level and assuming that 21 yum feature extends from 
18 to 22/rni we estimate the energy emitted in 21 ^m band 
as about 5.7 % of the total IR flux (251 L Q for A's from 5 to 
300 ^m assuming a distance to the source of 1 kpc). With the 
dotted line for wavelengths longer than 18 ^m we present the 
fit which was obtained using an opacity function with the EOF 
for only 21 /im component. It is clear that such fit is not able 
to explain IRAS photometry at 25 [im. Our recent ISO obser- 
vations show that this source is also a 30 /im emitter. In the 
forthcoming paper (Szczerba et al. 1999) we will discuss this 
finding in detail. 



Table 2. Model parameters for IRAS 04296+3429. Precise 
meaning of the symbols used can be found in Szczerba et 
al. (1997) 



parameter 


value 


T c ff 


6500 K 


log (L st ar [L Q ]) 


3.92 


d 


5.4 kpc 


Rout 


0.5 pc 


\r 

y cxp 


1 Vrrt c — 1 

iz Km S 


-R in (hot dust shell) 


6.4 10~ 4 pc 


T d [Rin(hotdust shell)] 


870 K 


pgas ^noi uusi snenj 


,-2.0 


Mpost-AGB 


4.010~ 7 M© yr" 1 


i?i n (main shell) 


7.06 10~ 3 pc 


T d [R in (main shell)] 


270 K 


p gas (main shell) 


r -2.6 

~ r 


jiVmin 
M AGB 


1.7O1O~ 5 M yr" 1 


iVmax 
M AGB 


2.19 lO^M© yr- 1 


Ct_ 


5A 


a+ 


0.25 ixm 


P 


3.5 


idyn 


575 yr 


M dust 


0.0071 M 



5. Discussion 

As it is shown in Table 1 , the metallicity for IRAS 04296 is 
significantly decreased relative to the solar value: the average 
abundance for the elements of the iron-group with respect to 
the Sun is [(Ti, V, Cr, Fe)/H] Q = -0.9 with the standard de- 
viation (7 = 0.2. 

Recently Decin et al. (1998), using high resolution spectra 
and model atmospheres method, calculated abundances of 14 
chemical elements in the IRAS 04296 atmosphere. Their re- 
sults are in qualitative agreement with these ones presented 
here, but there are some significant differences. Decin et 
al. (1998) calculated chemical composition of this object as- 
suming T cff = 7000K, log g= 1.0, £ t =4km/s, rather different 
from the model atmospheres parameters found in this work. 
It should be noted that we estimated the effective tempera- 
ture by two independent methods, and it is worth stressing that 
we have obtained consistent values of the effective tempera- 
ture: T ff = 6300 K from numerous Fe I, Fe II spectral lines and 
T ofi = 6500 K from modelling of the spectral energy distribu- 
tion of this source. The difference in effective temperature be- 
tween Decin et al. (1998) and our estimation (A T off = 700 K) 
is able to explain different metallicities estimated by Decin et 
al. (1998) and by us (Aloge(Fc) w 0.2). The same is true 
for the case of the rare-earth element abundances: large dif- 
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Table 3. Relative abundances of chemical elements for 
IRAS 04296+3429 in comparison to related PPNe. The data by 
Grevesse & Noels (1996) are adopted for solar abundances. 



IRAS 04296+3429 
[Fe/H] = -0.84 



IRAS 07134+1005" 
-1.00 



ROA 24" 
-1.77 



Element 



[X/Fe; 



o 



Li I 


> +0.23 






CI 


+0.84 


+1.08 


+0.67 


NI 


+0.83 


+1.03 


+1.02 


OI 


+0.19 


+0.63 


+1.01 


Nal 


+0.42 


+0.54 


+0.71 


Ms I 




+0.97 


+0.31 


Mali 


+1.34 




+0.09 


All 


+1.03 


+1.48 




Si I 


+0.58 


+0.95 


+0.80 


Sill 


+0.26 




+1.03 


SI 


+0.43 


+0.63 




Cal 


+0.19 


+0.45 


+0.60 


Sell 


+0.18 


-0.07 


-0.13 


Till 


-0.27 




+0.33 


VII 


+0.10 


-0.03 


+0.15 


CrII 


+0.11 




+0.65 


Cul 


+0.24 


+1.03 


-0.01 


Znl 


+0.08 






YII 


+1.20 


+1.70 


+0.37 


Zrll 


+0.62 






Ball 


+2.49 


+0.99 


+0.96 


La II 


+1.17 


+1.59 


+0.54 


Cell 


+0.82 




+1.60 


PrII 


+0.74 






Ndll 


+1.07 


+ 1.30 


+0.67 


EuII 


+0.34 


+1.06 


+0.25 



a -Klochkova (1995), 

b - Gonzalez & Wallerstein (1992). 



ferences, about 1 dex, in the values could be explained by dif- 
ferences in model atmosphere parameters. 

Let us consider now in more detail the peculiarities in the 
chemical composition of the object. For this purpose, in Table 3 
we present the logarithmic differences 

[X/Fe] = poge(X) - loge(Fe)]* - [loge(X) - loge(Fe)] 

between chemical compositions of different objects and the 
Sun (solar abundaces from Grevesse et al. (1996)): in the 
second column for IRAS 04296, in the third column for 
IRAS 07134+1005 (hereafter IRAS 07134 - asscociated with 
the peculiar F-type supergiant HD 56126) and in the fourth one 
for the star ROA 24. The objects are similar from point of view 
of their atmospheric parameters (T ff , log g) and relative chem- 
ical composition. It should be noted that the metal deficient su- 
pergiant ROA 24 (Fehrenbach's star) belongs to the globular 
cluster u Cen and could be considered as a typical halo object 
in the post-AGB evolution stage. 

The carbon overabundance [C/Fe] Q = +0.8 (revealed 
from intensities of 21 absorption lines with the stan- 



dard deviation a = 0.46) and the enhancement of nitro- 
gen [N/Fe] = +0.8 (from 4 lines, a = 0.10) suggest that 
IRAS 04296 underwent the third dredge-up episode. 

The oxygen content based on intensity of 3 weak lines near 
-A w 6155 A is determined with a small internal error. 

From the Fe-deficiency and CNO abundances (C/O > 2) 
_we can conclude that IRAS 04296 is a low mass object in ad- 
vanced stage of evolution. For an unevolved metal-deficient 
object (with [Fe/H] k, —0.9) the average value of [C/Fe] 
is only about -0.2 (Tomkin et al. 1995), the average value of 
[N/Fe] is w (Wheeler et al. 1989, Timmes et al. 1995) and 
the average value of [O/Fe] is « +0.5 (Wheeler et al. 1989, 
Timmes et al. 1995, Klochkova & Panchuk 1996b). The at- 
mospheres of the post-AGB stars IRAS 07134 and ROA 24 
are also overabundant in both carbon and nitrogen. Note how- 
ever, that for most of the PPNe candidates studied, strong rel- 
ative changes between elements of the CNO-group are ob- 
served (Luck et al. 1983; Lambert et al. 1988; Klochkova 1995; 
Zacs et al. 1995, 1996; van Winckel et al. 1996a, 1996b; van 
Winckel 1997). 

The abundances of some light metals (Na, Al, Mg, Si, 
Ca) are enhanced for all three stars. The average value for 
these elements is [X/Fe] = +0.6 for IRAS 04296; +0.9 for 
IRAS 07134 and +0.6 for ROA 24, with the standard devia- 
tions: a = 0.4, 0.4 and 0.36, respectively. 

We did not still include the K I abundance into our results, 
since we suspect that the equivalent width of its line near A 
7699 A could be significantly distorted due to circumstellar and 
interstellar components. 

The iron-group element zinc is the most important for de- 
termination of real (initial) value of the metallicity of a star 
"since, firstly, its abundance follows that of iron in a wide 
[Fe/H] interval (Sneden & Crocker 1988; Wheeler et al. 1989, 
Sneden et al. 1991) and, secondly, zinc having a low conden- 
sation temperature is not depleted by selective separation pro- 
cesses onto dust grains (Bond 1992). A close to solar abun- 
dance of Zn relative to iron ([Zn/Fc] = +0.1) permits us to 
conclude about the inefficiency of the selective separation pro- 
cesses in the IRAS 04296 envelope. This conclusion is based 
also on an absence of overdeficiency of light depleted elements 
(Ca, Sc). Besides, the relative abundance ([S/Fe] = +0.4 
with the standard deviation a = 0.21) of S, a chemical element 
which is not depleted by dust-gas separation, for IRAS 04296 
is close to the value for unevolved metal-deficient dwarfs 
(Francois 1987, Timmes et al. 1995). This futher confirms the 
lack of selective separation in the envelope of the object stud- 
ied. 

Individual abundances of the heavy s-process metals Y 
and Zr are determined with a relatively large error because 
of the small number of lines measured. However, the aver- 
age value [X/Fe] = +0.9 for Y and Zr is sufficiently reli- 
able. In addition, the abundance of heavy s-process element 
Ba ([Ba/H] = +2.5) derived from the equivalent width of 
strong lines could be altered by a systematic error due to the 
complexity of the outer regions of the stellar atmopshere as 
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discussed above. Nevertheless, we conclude that there is a Ba 
excess. 

The abundance of lanthanides (La, Ce, Pr, Nd) are 
strongly enhanced relative to iron for the objects from 
Table 3. For these heavy metals the average value is 
[X/Fe] Q = +1.0, +1.4, +0.9 for IRAS 04296, IRAS 07134 
and ROA 24, respectively, with the standard deviations 0.2 and 
0.6 for IRAS 04296 and ROA 24. Moreover, for all these ob- 
jects we see the overabundance of Eu which is predominantly 
produced by the r-process. 

Excess of s-process elements has been reliably found 
up to now in three objects investigated at the 6 m 
telescope: IRAS 04296+3429, IRAS 07134+1005 and 
IRAS 22272+5435. Besides, similar conclusions have ap- 
peared for another four PPN candidates (and for one object in 
common): HD 158616 (van Winckel et al. 1995); IRAS 19500- 
1709 = HD 187885 (van Winckel 1997); IRAS 05341+0852 
(Reddy et al. 1997); IRAS 22223+4327 and IRAS 04296+3429 
(Decin et al. 1998). In atmospheres of most PPN candi- 
dates overdeficiency (with respect to their metallicity) of 
heavy nuclei is generally observed (Klochkova 1995; van 
Winckel et al. 1996a, 1996b; Klochkova & Panchuk 1996a; 
van Winckel 1997), whose existence in the atmospheres of 
post-AGB low-mass supergiants has not yet found a clear 
explanation. 

In consequence, we can state that chemical abundances pat- 
tern for the source IRAS 04296 is related to its old galactic pop- 
ulation membership and dredge-up of matter enriched by the 
nucleosynthesis products. It may be part of the old disk popu- 
lation. 

As has been concluded already by Decin et al. (1998) all 
the post-AGB candidates mentioned above (only these, up to 
now, show an s-process element enhancement!) belong to the 
small group of PPNe (Kwok et al. 1989; Kwok et al. 1995) 
which have in their IR spectrum an unidentified emission band 
at about 21 /mi. This feature is neither found in the spectra 
of their predecessors, AGB stars, nor in the spectra of PNe. 
Note, once more, that the search by means of the ISO for the 
new 21 fjm emitters among candidates selected by Henning 
et al. (1996) failed (Henning, private communication). As has 
been stated in the papers by Kwok et al. (1989, 1995), the ob- 
jects whose spectra contain the 21 fira band are carbon-rich 
stars. Our investigations based on the spectra from the 6 m tele- 
scope, for IRAS 07134 (Klochkova 1995), IRAS 22272+5435 
(Zacs et al. 1995) and IRAS 04296 (Klochkova et al. 1997b), 
confirmed that C/O > 1 for all of them. In this context, the 
conclusion that the carrier of the 21 /mi band is related to C is 
natural. For example, Buss et al. (1990) have supposed that this 
feature may be caused by polycyclic aromatic hydrocarbons. 
On the other hand, Goebel ( 1 993) has identified the 2 1 /im band 
with the vibrational band of the SiS 2 molecule, the presence of 
which is consistent with the temperature in the envelope. 

Taking into account the available results on chem- 
ical composition for subclass of PPNe with the 
21 /mi feature: IRAS 07134+1005 (Parthasarathy et 
al. 1992, Klochkova 1995), IRAS 22272+5435 (Zacs et 



al. 1995), IRAS 19500-1709 (van Winckel et al. 1996a), 
IRAS 05341+0852 (Reddy et al. 1997), IRAS 22223+4327 
(Decin et al. 1998), and IRAS 04296 (Decin et al. 1998; 
Klochkova et al. 1997b; this paper) we see that the carbon-rich 
atmospheres of these objects are also enriched by s-process 
elements. It is evident that there is a strong correlation 
between presence of the 21 /im feature, C2, C3 molecu- 
lar bands, and excess of the s-process elements. Decin et 
al. (1998) were the first who pointed out this relationship. 
What is even more important, an excess of s-process ele- 
ments was not found for a number of IRAS sources with 
altered CNO-content but without the 21 /mi feature (some of 
which are oxygen-rich stars rather than carbon-rich stars): 
IRAS 06338+5333 (Luck & Bond 1984; Bond & Luck 1987), 
IRAS 07331+0021 (Luck & Bond 1989; Klochkova & 
Panchuk 1996a), IRAS 09276+4454 (Klochkova & Mishen- 
ina 1998), IRAS 12175-5338 (van Winckel 1997), IRAS 
12538-261 1 (Luck et al. 1983; Klochkova & Panchuk 1988b; 
Giridhar et al. 1997), IRAS 15039-4806 (van Winckel et 
al. 1996b), IRAS 17436+5003 (Klochkova & Panchuk 1988a; 
Luck et al. 1990; Klochkova 1998), IRAS 18095+2704 
(Klochkova 1995), and IRAS 19114+0002 (Zacs et al. 1996, 
Klochkova 1998). Therefore, it seems that carrier of 21 /im 
feature is strongly related to the whole chemical composition 
pattern typical for the third dredge-up (excess of s-process 
elements), and not only to the C-richness of the photosphere. 

That 21 /mi feature is not observed around AGB-stars 
showing s-process elements could be explained by the physi- 
cal conditions which are inappropriate for the excitation of this 
band, while its non-presence in planetary nebulae may be a 
result of carrier destruction by the highly energetic photons. 

6. Conclusions 

We conlude that IRAS 04296+3429 is a PPN candidate with 
a chemical composition which coincides with theoretical pre- 
dictions for the post-AGB objects: very large excess of carbon 
and nitrogen are revealed. Moreover, the real excess (relative to 
iron) of heavy metals Y, Zr, Ba, La, Ce, Pr, Nd synthesized by 
the neutronization process indicates an effective third dredge- 
up and further confirms IRAS 04296 to be in the advanced post- 
AGB evolution stage. 

The emission of C2 molecular lines discovered in the spec- 
tra of IRAS 04296 and its similarity to the emission of the 
Hale-Bopp comet allow us to suggest that in both cases the 
same mechanism (the resonance fluorescence) is responsible 
for the observed features. 

Several strong absorptional features whose positions coin- 
cide with known diffuse interstellar bands (DIBs) are found in 
the spectrum of IRAS 04296. 

In addition, from the SED modelling of the spectral en- 
ergy distribution we showed that 25 fxm flux cannont be ex- 
plained without assumption that the 30 /im emission feature is 
present in this source. Our recent ISO observation has detected 
the 30 /im band in this source. 
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